Field-induced quasi-periodic coherence effects in certain small magnets. 
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Small ferromagnets and anti-ferromagnets with an easy-plane anisotropy have a ground to first ex- 
cited state (tunnel) splitting which is quasi-periodic in the magnitude of a field applied perpendicular 
to a principal anisotropy axis. The associated oscillations in thermodynamic quantities might be 
used to prove the existence of a coherent ground state even when the tunnel splitting itself cannot 
be directly detected. 
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The study of small magnets tests our understanding of 
the transition between quantum and classical physics. In 
the classical regime the order parameter, i.e., the mag- 
netization, of a ferromagnet points is one of a number of 
equivalent directions which correspond to a minimum of 
the macroscopic anisotropy energy. On the other hand 
in the quantum world of small magnets the ground state 
comprises a coherent superposition of all of these equiv- 
alent equilibrium directions. While there is now fairly 
strong evidence JO of relaxation via quantum tunneling 
in Mni2 the observation of a coherent ground state re- 
mains an elusive goal. The only experimental evidence, 
for the anti- ferromagnet ferritin 0, remains highly con- 
troversial. 

The principal purpose of this Letter is to show, for 
small ferromagnets and anti-ferromagnets with an easy- 
plane anisotropy, there exist rather dramatic level cross- 
ing effects. Unlike similar easy-axis magnets which have 
tunneling effects very sharply peaked at level crossings, 
the tunnel splitting oscillates as a function of the mag- 
net field. The period, e.g., AH — 2H» for a ferromag- 
net, is determined by an easy-to-measure bulk anisotropy 
field Hh . The observation of this oscillatory signature 
can provide an unrefutable proof of coherence, and since 
M = —dF/dH where F is the free energy, there is neces- 
sarily an oscillatory component in the magnetization, and 
other thermodynamic quantities, which might be readily 
detected even if the splitting cannot. 

The ferromagnet is modeled by a single large spin sub- 
ject to the external and anisotropy fields fjl, i.e., the 
Hamiltonian Ti — gfisS ■ H + K\\S Z 2 + K±S X 2 where 
without loss of generality it is assumed that > \K±\. 
For an easy plane magnet K\\ > 0. The energy parame- 
ters Kn, K± and, for the anti- ferromagnet, J all scale as 
S . The equivalent physical quantities are Hu = Kn/ S, 
H± = Kj_/S and H e = J/S. This is important in ratios 
such as (h/K l{ ) = (Sh/H\\). 

The problem is formulated |||J in terms of auxiliary 
particles. Consider a ferromagnet. A basis \S Z >= \n > 
is chosen. Then an auxiliary particle, a fermion f n , is 
associated with each state via the mapping \n >— > > 
where | > is a non-physical vacuum without any auxiliary 



particles. Defined is a bi-quadratic version of an operator 
6 via: 6 — > J2nn>fn < n \0\n' > f n ,. The constraint 
Q = J2 n = J2 n fnfn = 1- It has been shown [| that 
such schemes preserve all operator multiplication rules 
including commutation rules. The replacement rule is 
applied to the Hamiltonian Ti. to yield, taking gfJ-sH = 
—hi: 

\K ± [(M^) 2 + (M«- 1 ) 2 ])/t/ n 
+ \K ± £„ M^MZttUUzU + H.c), 



(1) 



where the M™ +1 = [S(S + 1) - n(n + 1)] 1/2 are the ma- 
trix elements of S' ± . This is two tight binding models of 
spinless fermions f^. The constraint Q = 1 implies this 
is a single particle problem. 

The two "chain" structure reflects the fact that the 
"hopping" term in Eqn. ([j]) couples "sites" with indices 
which differ by two. This structure implies immediately 
a spin-parity effect found by Loss et al. ||, and von Delft 
and Henley |Q . For the case of integer spin, Fig. |l|, the 
two chains comprise the sites n = — S, — (S — 2),— (S — 
3)...(S-3),(S-2),S&ndn = -(5-1), -(5-3), -(S- 
5), ... (S — 5), (S — 3), (S — 1) which are both symmetric 
relative to n = 0. On the other hand for half-integer 
spin, Fig. ||, the chains are n = —S,—(S — 2),— (S — 
3) ... (S - 5), (S - 3), (S - 1) and n = -(S - 1), -(S - 
3), —(5 — 5), ... (5 — 3), (5 — 2), 5, which are equivalent to 
each other through the map n — > — n but which lack the 
symmetry about n = 0. Because of the equivalence of 
the two chains there must always be a double (Kramers') 
degenerate ground state, without a tunnel splitting, for 
half-integer spin, while because of symmetry about n — 
0, there can be tunnel split pairs for integer spin. 
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n= -5-3-1 13 5 

FIG. 1. The two chains for integer spin each have symme- 
try about n = 
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n=-13/2 -9/2 -5/2 -1/2 3/2 7/2 11/2 

n= -11/2 -7/2 -3/2 1/2 5/2 9/2 13/2 

FIG. 2. The chains for half-integer spin are equivalent 
though the map n «-» — n 

There is a key /ixerf pom£ [p| when K± — > 0. All off- 
diagonal matrix elements are zero, the diagonal energies 
are K^n 2 — hn, and for an integer spin easy-plane fer- 
romagnet, i.e., with K\\ > 0, and in zero field (h = 0), 
the ground state lies on the site n = 0, and therefore the 
even site chain. The first excited states at an energy Kn 
are located on the sites n = ±1 and corresponding to the 
odd site chain. Level crossings occur whenever 

h = (2n + l)Kf n = 0,1,2,3..., (2) 

which has the period 2Ku . With this period the splitting 
between the ground state and the first excited state un- 
dergoes sawtooth oscillations with amplitude ~ Ku and 
with each level crossing the ground state passes from one 
chain to the other. Given that as a function of K± there 
are no singularities in the characteristic determinant, this 
fixed point will govern the behavior for all values of this 
parameter. That this is the case for a ferromagnet is 
verified by both the analytic and numerical calculations 
described below. 

The vicinity of the fixed point with 

S 2 K ± < K\\ (3) 

defines the small particle limit. The corrections to the 
ground state \n — > are perturbative, however the ex- 
cited states split. These are (l/v / 2)[jl > ±| — 1 >] with 
energies K«±K±. The smallest h — splitting K\\ — \K±\ 
determines the amplitude of the near saw-tooth oscilla- 
tions. This splitting has nothing to do with tunneling. 
However the ground state is clearly "coherent" in the 
sense < S > = 0, i.e., the order parameter is not local- 
ized in a particular direction in the x — y-plane. And if, 
e.g., the initial state is (1/V3)[|0 > +|1 > +| - 1 >], 
< S > has a magnitude of ~ S and points in the posi- 
tive x-direction. With a frequency io\ = Hu , it oscillates 
rapidly between the +x and — ^-directions. While with 
the frequency L02 = H±_ these oscillations precess in the 
perpendicular plane so that one quarter of a long period 
later the rapid oscillations are between the ±y-directions. 
Such oscillatory behavior is characteristic of coherence. 
Quite generally Schrodinger's equation 

(e — (K^n 2 — nh))a n 

= ±iT x M^ 1 Af£+?[(a tl+ 2 - a n ) + (a„_ 2 - a„)], (4) 

involves finite differences, where the wave-function >F = 
J2n(~ )" a «/nl >• When the inequality (||) is re- 
versed the particle is large and a continuum approx- 
imation to this Schrodinger's equation is appropriate. 



Within the definition of a Fourier transform, a n = 
(1/V2tt) /J™ dpf{p)e~ tp< - n - d \ the function f(p) is de- 
fined for the even site chain and h = 0. For finite h 
the center of the real space potential is displaced by 
d = h/2K\\ and the odd site chain is accommodated via 
a d = 1 + {h/2K\\). Assuming K± < 0, the continuum 
approximation to Eqn. (Q) is: 

{e -^ +K W 2)s{p) = -^t^- 1 ^ 

(5) 

which is Mathieu's equation ||. The potential is pe- 
riodic with period Ap = 7r and the low lying so- 
lutions form bands. For a given energy, a solution 
might be characterized by a wave- vector k and Floquet's 
(Bloch's) theorem S implies that solutions are of the 
form f k (p) = e lkp u k {p) where u k (p) = u k (p + n). That 
a n be finite implies k — d (to within trivial translations). 
Within a given well of the periodic potential, the low 
lying states have wave-functions which to a good ap- 
proximation are those of a harmonic oscillator. Taking 
K x < 0, around p = 0, f(p) = (1/ 1 yJJ^)e~ p ' 'W with 
P 2 = Sy/Ku/Kj_, and the nominal ground state energy is 
(h 2 /4K\\) + (Lu Q /2)- lu = 2Sy/K\\K±. Because of tunnel- 
ing between wells this ground state level becomes a band, 
i.e., the energies e k = (/i 2 /4if||) + (wq/2) + (to/2) coswk 

where the width || w = 8y /, 2/i:ujoSy 2 e~ Sf and where 
the action S f = 2S y J{K ± /K\ ] ). 

The result for the "tunnel splitting", i.e., the difference 
in energy between the ground states for the even and odd 
chains is 

^ = 4^04/^^003^); 

S f = 2Sy/(HJH\\). (6) 
The numerical work shown in Fig . confirms this result 
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FIG. 3. Each plot corresponds to K\\ = 1. For (a) 
K± = .0003, S — 40, and the saw-tooth tunnel spilting 
~ if II has been divided by 10 4 . Both (b), with integer spin 
S = 40, and and (c), with the half-integer value S = 81/2, 
have K± = .03 and the form is sinusoidal. Notice the ir/2 
phase shift. 
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that large particles have a tunnel splitting which oscil- 
lates as cos(S , 7r/i/2_ff||). This with the equivalents for an 
anti-ferromagnet, are the principal results of this Letter. 

A perpendicular magnetic field, i.e., one lying in the 
x — y-plane, couples the two chains. Still, for a larger 
particle and Kj_ < 0, the wave-function, a n , is ex- 
tended, does not alternate in sign between sites, and is 
well localized near the center of a chain. The operators 
5 ± , to a good approximation, cause translations by one 
lattice site, i.e., they convert the ground state on one 
chain into that on the other. The matrix elements of 
h x S x — (h x /2)[S + + S~] are approximately h x S/2 while 
those of hySy are negligible. Assuming that h is smaller 
than the splitting ~ luq between the the ground doublet 
and other excited states, the splitting of the doublet is: 



SE h = ^(SE) 2 + (Sh x y 



(7) 



Evidently, the field scale it determined by the tunnel 
splitting SE/S and implies a rather careful alignment is 
required in order to observe oscillations. If K± > 0, a n 
alternates in sign between sites and the role of the x and 
y axes are interchanged. 

Again for a large particle (and K± < 0), since S x has 
large, and S y very small, matrix elements within the low 
lying doublet, it follows that if the initial condition puts 
the system in an approximately equally weighted linear 
combination of these two states, the magnetization will 
oscillate in the cc-direction at a frequency correspond- 
ing to the tunnel splitting, but the component in the 
y-direction is negligible. Thus, in passing from a small 
to a large particle the magnetization becomes localized 
in the x-direction. Again for K± > the role of the x 
and y axes are interchanged. 

Turning to spin-parity effects, consider a small parti- 
cle with half-integer spin. The fixed point now has level 
crossings occurring whenever h = 2nKu n = 0, 1, 2, 3 . . . 
which reflects a shift of tt/2 relative to the integer spin 
case. Kramers' degeneracy corresponds to the level cross- 
ing at h = 0, but apart from the tt/2 phase shift a half 
integer small particle behaves like its integer spin equiv- 
alent. 

Similarly for a large particle and half-integer spin, all 
is the same as for the integer spin case except for the tt/2 
phase shift. In particular in the tunnel splitting, Eqn. 
(g), the cosine is replaced by a sine. 

An easy- axis ferromagnet [^) is quite different from its 
easy-plane equivalent. In particular there is no oscilla- 
tory field dependence. 

For larger fields an easy-plane anti-ferromagnet is 
much like a similar ferromagnet and there is again a 
quasi-periodic modulation of the tunnel splitting. The 
effective Hamiltonian [[TJJ H = g^ B S ■ H + JSa ■ Sb + 

K\\ (Sa,z) 2 + (Sb.z) 2 contains the external field, ex- 
change and a suitable anisotropy energy. The exchange 
field J couples the sub-lattices 11 A" and "_B" and it is 



assumed that Sa = Sb = S unless stated otherwise. 
The auxiliary particles a n m create single particle states 
\n, m >= a\ lm \> which map to the \Sf = n + m, = 

—n >. The constraint is Q m — Y^n. a nm a nm = 1j anc ^ 
the auxiliary particle Hamiltonian 



E 



\ a n+l. : 



+ H.c.) (8) 



where, the diagonal energies are e n = —J(n + m)n — 
Kii[(n + m) 2 + n 2 ], and the hopping matrix elements 

C +1 = { J l 2 ) M n+m +lM n +l - Hopping couples neighbor- 
ing sites, and there is a distinct chain for every value of 
the quantum number m = S z . 

The chains are roughly equivalent to the sites for the 
ferromagnet. The chains with m values which differ by 2 
are coupled by the inter-chain hopping terms which arises 
when the term K±\Sa, x 2 + Sb, x 2 ] is added. The struc- 
ture comprises two distinct two dimensional networks, 
the equivalent of the chains for a ferromagnet, and as 
the magnetic field changes the ground state can pass form 
one network to the other causing the tunnel splitting to 
oscillate in a quasi-periodic manner. Here it is not as- 
sumed that K± < Kn . The larger and smaller of the two 
will be denoted K and k respectively. 

The period of the oscillations is determined by the con- 
dition for level crossings. The details are different for 
small and large particles. A small particle is defined by 



J > S 2 Ku 



(9) 



For a given chain, with index m (= S z ), the ex- 
change energy JSa • Sb dominates the anisotropy en- 
ergy K» and trivially the ground state has an energy 0] 
(J/2)|m|(|m| + 1) — mh. The fixed point level crossings 
occur when 



h=J(n+l); n= 1,2,3. 



(10) 



i.e., the first possible crossing, and a zero of the tunnel 
splitting, occurs when h = 2J and otherwise has period 
Ah = J. Notice there is a phase shift of tt/2 relative 
to level crossings for a similar ferromagnet, Eqn(^), and 
that J/2 has replaced the anisotropy parameter K\\ . The 
magnitude of the tunnel splitting depends upon another 
small to large cross-over. If S 2 K± < J the particle is 
twice times small and the periodic saw-tooth like split- 
ting is ~ J = H e /S. Here the anisotropy energy is irrel- 
evant and similar periodic behavior will occur for a twice 
small easy axis anti-ferromagnet. 

When this latter inequality is reversed the particle is 
small-large and it is implied that K± = K = H a /S, 
Kh — k. There are two regimes. For the small field 
regime, h < uq 



2S\/KJ = 2\/H a H e , the axis of quan- 
tization is rotated to lie along the x-direction, which in- 
terchanges the role of, K± = k {= for simplicity) and 
K\\ = K and the field term is hS x . As a result, the 
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inverse of inequality (Q) is satisfied, the wave-functions 
on the chains are those of an harmonic oscillator, and 
the chain ground state energies are Q (ojq/2) + Jm 2 /2. 
The very small field limit applies when h < J = H e /S 
whence the coupling between the chains is a perturba- 
tion. The field couples each chain m to the first excited 
state on two adjacent chains mil, with a matrix el- 
ement h(3/2*J2, P 2 = 2S y / K/J. The diagonal energy 
correction is — 2(/i/3/2v / 2) 2 / CJ o- An off-diagonal term 
th = — (hf3/2\/2) 2 /luq couples the m = ±1 chains via 
an excited state on the m — chain. This splits the pair, 
with the lower energy being (J/2) — 3(h(3) 2 /8uJo- For very 
small fields the splitting is therefore 



(hp) 2 
8luq 



J 
2 



8T 



SH, 



(ID 

8 H R v ' 



For larger fields, but still in the small field regime, the 
off-diagonal coupling th exceeds the difference in energy 
between adjacent sites and a continuum approximation 
is appropriate. The problem reduces |l2| to that for a 
ferromagnet and the result is Eqn. (|6|) with h = 0, 



#11 -> H e /2 and H± 



relevant experimental evidence known to the author is 
for the anti-ferromagnet Feio 10 an d would appear to 
belong to the small-small limit. The series of small fer- 
romagnets based on Mn, all seem to have an easy axis. 
However many other systems exist for which the nature 
of the anisotropy have never been determined. 
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-h /2H e . The zero point en- 
ergy for motion perpendicular to the chains is ~ h while 
that associated with motion along the chains is vq. The 
equality of these two energies indicates the point at which 
the appropriate axis of quantization changes and implies 
the small to large field inequality quoted above. 

In the large field regime, h > luq, reverting to an 
axis of quantization defined by the field direction, the 
ground state energy for the chain labeled m is again 
(J/2)\m\(\m\ + l)-mh= (J/2)[m - (h/J) + (1/2)] 2 - 
(h 2 /2J) — (J/8) and represents an harmonic potential, 
centered at {h/J)— (1/2), for the motion perpendicular to 
the chains. The coupling between chains is S 2 K±/2 and 
Eqn. (||) again applies with ify — > H e /2, H± — ► 2H±, 
and with the cosine replaced by a sine p2[ . The same 
result is valid for small integer values of j = \Sa — Sb\ 
the excess, or net, spin. For small half-integer j values 
after the substitution H\\ — > H e /2, Eqn. (|J) is valid, as 
written. 

For a larger magnet when inequality (|^) is reversed, the 
chain ground state energy is again (luq/2) + (l/2)m 2 J — 
mh, identical, apart from a constant, to the site en- 
ergy for the ferromagnetic problem with the replace- 
ment if 1 1 — > J/2. The level crossings correspond to 
h = J(2n + l)/2, i.e., have period J. Also as for the 
previous cases, the magnitude of the tunnel splitting de- 
pends upon another small to large cross-over. Now if 
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half-integer j. 

Suitable single crystals of molecular magnets would 
ideal for studying the predictions made here. The only 
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